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Abstract 

In this article we derive longitudinal impedance and wake function for an undula- 
tor setup with arbitrary undulator parameter, taking into account a finite transverse 
size of the electron bunch. Earlier studies considered a line density-distribution of 
electrons instead. We focus our attention on the long-wavelength asymptote (com- 
pared with resonance wavelength), at large distance of the electron bunch from 
the undulator entrance compared to the overtaking length, and for large vacuum- 
chamber size compared to the typical transverse size of the field. These restrictions 
define a parameter region of interest for practical applications. We calculate a closed 
expression for impedance and wake function that may be evaluated numerically in 
the most general case. Such expression allows us to derive an analytical solution for 
a Gaussian transverse and longitudinal bunch shape. Finally, we study the feasibil- 
ity of current-enhanced SASE schemes (ESASE) recently proposed for LCLS, that 
fall well-within our approximations. Numerical estimations presented in this paper 
indicate that impedance-induced energy spread is sufficient to seriously degrade 
the FEL performance. Our conclusion is in contrast with results in literature, where 
wake calculations for the LCLS case are given in free-space, as if the presence of 
the undulator were negligible. 
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1 Introduction 



Self-Amplified Spontaneous Emission Free Electron Lasers (SASE FELs) 
IHLI2II3I] are nowadays considered as a unique tool for production of intense, 
polarized, short-pulse radiation tunable throughout the VUV and X-ray 
wavelength range, with peak and average brilliance exceeding both modern 
Synchrotron Radiation and Laser Plasma sources by many order of magni- 
tudes |0,|5i@,0r&@]. Successful operation of SASE FELs requires high quality 
(low emittance and low energy spread), intense electron beams. One of the 
trends for SASE FELs is production of ultra-short radiation pulses. These 
can be obtained by exploiting electron bunches with an ultra-short charge 
concentration (spike). At the FLASL-Q facility at DESY, Hamburg 
electron bunches with sharp spikes have been produced in the framework 
of a nonlinear bunch-compression scheme. Experimental and 
theoretical |l(i H] studies of FLASH operation have shown that properties 



of ultra-short pulses are significantly influenced by collective effects, the 
most important of them being space-charge effects. Space-charge plays an 
important role in the beam-formation system, in the drift space and also in 
a long undulator. Collective effects might be crucially important for X-ray 
SASE FELs (XFELs) as well. 

This article presents a description of longitudinal wake fields in XFELs. In 
particular, our study is of importance in connection with novel schemes 



of radiation production, like Enhanced SASE schemes (ESASE) |12l, 1131, 114(1 . 
ESASE proposals rely on two steps. First, the electron beam is modulated in 
energy by interacting with a GW-level optical laser in a modulator wiggler 
placed in the accelerator section. Second, a dispersive section transforms 
the energy modulation into density modulation, eventually leading to a 
subfemtosecond-long spike in the beam current before the entrance in the 
FEL undulator. The peak current of this spike can reach tens of kA without 
emittance worsening, because only a small charge is concentrated in the 
high-current region. As the electron beam undergoes the SASE process, the 
enhanced current part should saturate faster than the rest of the bunch. 
Alternatively, the x-ray wavelength may be reduced, for a fixed undulator 
length. Moreover ultra-short pulses (in the attosecond range) are produced 
as a result of the presence of the short lasing spike. Faster saturation of 
emission from the enhanced-current spike also suggests that ESASE schemes 
may be used to obtain saturation even in situations when beam parameters 
deteriorate with respect to design values. 

A detailed study of longitudinal wake fields arising after the dispersive 
section, in particular dominant space-charge wake fields is due in order to 



Free-Electron Laser in Hamburg. 
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assess the magnitude of detrimental effects on the FEL process. It is impor- 
tant to note that the undulator parameter K for XFEL setups obeys K 2 » 1. 
As a result, the average longitudinal Lorentz factor y z = yj Vl + K 2 /2 is 
such that y\ y 2 , y being the Lorentz factor of the beam. Based on y\ «: y 2 , 
we will demonstrate that the presence of the undulator strongly influences 



the space-charge wake. In contrast to this, in ||13lll4{1 , wake calculations for 
the LCLS case are given in free-space, as if the presence of the undulator 
were negligible. Authors of references 1 13 r lU] incorrectly conclude that the 



FEL process is basically unaffected by space-charge wakes. 

This paper is devoted to the calculation of impedance and longitudinal 
wake field in XFELs, with particular attention to the LCLS case, for which 
ESASE schemes have been first proposed. This means that we will restrict 
our attention to a very specific region of parameters, discussed in the next 
Section 2. First, the longitudinal size of the beam is much larger than the 
FEL wavelength. Second, electrons are assumed to have travelled into the 
undulator for a distance longer than the overtaking length. Third, effects 
of metallic surroundings can be neglected. When the electron-beam size is 
larger than the radiation diffraction size calculated from a single undulator 
periocGl major simplifications arise. In fact, radiation from the undulator 
is drastically suppressed and calculations of impedance and wake function 
can be performed considering a non-radiating beam, and thus accounting 
for space-charge interactions only. Then, space-charge impedance and wake 
function is found to reproduce the free-space case. Only, the Lorentz factor 
y must be consistently substituted with the average longitudinal Lorentz 
factor y z . In Section 3 we derive the electric field that will be used to calculate 
impedance and wake. In Section 4 we introduce concepts of impedance and 
wake field. Fields are calculated in Section 3, while impedances and wakes 
are respectively dealt with in Section 5 and Section 6. Then, in Section 7, 
we apply our theory to the ESASE setup referring to the LCLS facility. We 
calculate the energy chirp associated with wakes inside the undulator and 
between dispersive section and undulator. Subsequently, the magnitude of 
their effect is estimated by calculating the linear energy chirp parameter 
115, 3]. We find that the gain of the FEL process is sensibly reduced, and 



that longitudinal wake fields constitute a reason of concern regarding the 
practical realization of ESASE schemes. Conclusions end our treatment in 
Section 8. 



2 Of order ^AA W , A being the reduced wavelength of coherent radiation and K u 
the undulator period. 
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2 Parameter-space of the problem 



As has been said in the Introduction, results of this paper can be applied to 
calculate effects of wake fields in planaiOJundulators under a specific choice 
of parameters corresponding to an XFEL system. Quantities of interest are 
defined once the bunch and the undulator system are specified. The bunch 
is characterized by an rms length a Z/ a transverse rms dimension o± and 
Lorentz factor y. Moreover, we define the undulator period \ w , the vacuum 
chamber transverse dimension a and the undulator parameter K, where 
K = A w eH w l(2nm e c 2 ), (—e) being the negative electron charge, H w the peak 
undulator magnetic field on-axis, and m e the rest mass of the electron. Finally, 
L s is the saturation length of the FEL process. 

The bunch length o z corresponds to the reduced wavelength of the coherent 
field generated by the bunch, a z A = A/(2n). This (reduced) wavelength 
is much longer than the reduced resonant wavelength A,- - A w / (2y 2 ), where 
Ao, = A w /(2n), and y z = yj Vl + K 2 /2 is the already defined average longi- 
tudinal Lorentz factor. This means A » A,.. 

The overtaking length is defined by the quantity 2y\ A. When the bunch has 
travelled inside the undulator for more than 2y\ A a steady state is reached, 
and asymptotic expressions for the wake fields can be given. In the present 
study we will work with such asymptotic expressions only. This means that 
the saturation length of the FEL process, L s , must be much longer than the 
overtaking length, i.e. L s » 2y\K. 

Also, in this paper we will neglect the presence of the vacuum chamber. This 
is possible when the vacuum chamber dimension is much larger than y z A, 
i.e. a » y 2 A, a typical transverse dimension associated with the coherent 
field, that is verified for ESASE XFEL setups. 

Summing up, we will work under the following constraints: 



A » A,. , 
L s » 2y\K, 

a » y z K . (1) 

Based on conditions in (1), we will develop a theory of wake fields from 
undulators in XFELs. In particular, the first assumption greatly simplifies 



Although we presented final expressions of our theory in the case of a planar 
undulator, there are no specific effects related with the choice of a planar undulator. 
Our work may be straightforwardly extended to the case of a helical undulator as 
well. 
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our consideration allowing for a long-wavelength asymptotic treatment. 
Under the second and the third assumption we will be able to present an 
expression for the impedance in terms of a double convolution involving 
the charge density distribution and Bessel functions. Similarly, an analytical 
expression for the wake function could be given. However, it will not be 
necessary to explicitly calculate this expression. In fact, when discussing 
practical applications, we will work in the asymptotic case 



o\ » KK W . (2) 

Extra-condition (2) greatly simplifies the treatment of wake fields. Our re- 
sults can be directly applied to realistic situations as the ESASE scheme 
analyzed in Section 7, where we will refer, explicitly, to the LCLS case. 



3 Field calculation 



Calculation of longitudinal wake field and impedance from an FEL undula- 
tor is subject to the characterization of the electric field generated at a given 
position (that is the position of a test electron) by the entire bunch. 

We perform an analysis in terms of harmonics, i.e. we consider sinusoidal 

dependence of the electric field of the kind E = E(f. co) exp[-icoi\ + C.C., 
the symbol "C.C." indicating complex conjugatioro Here t is the time, 
co = 2nc/A is the frequency, with c the speed of light in vacuum. The complex 

amplitude E(r,co) can actually be considered as the representation of the 
electric field in the space-frequency domain, and it will be referred to as 
"the field". 

We assume that particles proceed along an undulator, under the constraints 
discussed in Section 2. 



3.1 Transverse field 



The transverse field E ± can be treated in terms of Paraxial Maxwell's equa 



tions in the space-frequency domain (see e.g. [13, LL8Q). From the paraxial 



4 In the following, for simplicity, we will consider co > 0. Expressions for the field at 
negative values of co can be obtained based on the property E(-co) = E (co) starting 
from explicit expressions for E at co > 0. 
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approximation follows that the electric field envelope E ± = E ± exp [-icoz/c] 
does not vary much along z on the scale of the reduced wavelength K = 
A/(2n). As a result, the following field equation holds: 



E ± (z,r ±> co) 



g(z,f ±/ co) , 



(3) 



where the differential operator D is defined by 



»-(v^). (4) 

V ± 2 being the Laplacian operator over transverse cartesian coordinates. 
Eq. (3) is Maxwell's equation in paraxial approximation. The source-term 
vector g(z, r ± ) is specified by the trajectory of the source electrons, and can be 
written in terms of the Fourier transform of the transverse current density, 

j ± (z, r±, (o), and of the charge density, p(z, r ± , co), as 



g = -4tt exp |-— j - V ± p j 



(5) 



j ± and p are regarded as given data. In this paper we will treat j ± and p as 
macroscopic quantities, without investigating individual electron contribu- 
tions. We consider transverse and longitudinal distribution densities of the 
current constant through the undulator. In the time domain, we may write 
the charge density p(r, t) and the current density j(r, t) as 



1 P x(rx-rUz))/(f-^ 



V z(z) 



(6) 



and 



H?> f ) = — -T^v^p^ - r' 01 _{z))f t 



s (z) N 



V z(z) 



(7) 



The quantity p ± has the meaning of transverse electron beam distribution, 
while / is the longitudinal charge density distribution. r' 0± (z), s (z) and v 
pertain a reference electron with Lorentz factor y that is injected on axis 
with no deflection and is guided by the undulator field only. Such electron 
follows a trajectory specified by f OJ _(z) = r' ox e x + r' oy e y , e x and e y being the unit 
vectors in the horizontal and vertical directions respectively, with 
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C( z ) = -j^ cos(k w z) = r w cos(k w z) , r' oy (z) = , (8) 

where we defined the transverse amplitude of oscillations r w = K/ (yk w ). The 
corresponding velocity is indicated with v 0± _{z) = v ox e x + v oy e y : 



Kc 

v ox (z) = sin(k w z) , v oy (z) = . (9) 

y 

Finally, s (z) is the curvilinear abscissa measured along the trajectory of the 
reference particle. 

Note that, according to Eq. (6) and Eq. (7), / = v p. In fact, for each particle 
in the beam byly <sc 1. Therefore we can neglect differences between the 
average transverse velocity of electrons (v) and v . 

In the space-frequency domain, Eq. (6) and Eq. (7) transform to: 

p(r ± , z, w) = p (r ± - rV(z)) f(co) exp [icos (z)/v ] (10) 

and 

(z))/M 

exp [icos (z)/v ] , (11) 
where, for simplicity, we introduced the symbol 



P°fo = ^TmP^ ■ (12) 

Voz\Z) 

-± 

It should be remarked that p and ; = pv satisfy the continuity equation. In 

-> -± 

other words, one can find V • / = icop. 
We note that for a generic motion one has 

/ s (z 2 ) - S (Zi) z 2 -z 1 \ f CO 

"I — » — ) = J (13) 



Also, 



z 

/£l> £l> wK 2 . x CO 
— dz = tZ - — — — sm(2k w z) -z , (14) 
2cyl{z) 2cyl 8y 2 k w c v ' 2cf z V ' 

o 
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where the average longitudinal Lorentz factor y z is defined as 

?z = 7 ■ (15) 



The approximate equality in Eq. (14) follows from the fact that we are 
interested in wavelengths to <3C co r , where the fundamental co r = 7k w cy\ 
is fixed imposing resonance condition between electric field and reference 
particle. The term in sm(2k w z') is of order co/co r , that is our accuracy, and can 
be neglected everywhere. 

With the help of Eq. (10) and Eq. (11), Eq. (5) can be presented as (see also 
£3): 



p (j±-^V(z))/M . 

(16) 

We find an exact solution of Eq. (4) without any other assumption about 
the parameters of the problem. A Green's function for Eq. (4), namely the 
solution corresponding to the unit point source can be written as (see ilZy): 



g = -An exp 



z 

I 



dz- 



a) 



2f z c 



-tV 0± (z) - Vj 
L c z 



G(z-z';r:-ri) = - i ^- 7y exp < 

assuming z - z' > 0. When z - z' < the paraxial approximation does not 
hold, and the paraxial wave equation Eq. (3) should be substituted, in the 
space-frequency domain, by a more general Helmholtz equation. Yet, the 
radiation formation length for z - z' < is very short with respect to the case 
z—z' > 0, i.e. we can neglect contributions from sources located at z - z' < 0. 

Thus, after integration by parts, we obtain the solution 



ico- 



r, -r. 



2c(z - z') 



(17) 



2 



too 



r± - r' ± Y 
2c(z - z') 



I' 



+ i I dz 



CO 



'2cy 2 z (z) 



1 /_, x _ lv x (z') rj_-^j_\ 



(18) 



Eq. (18) describes the field at any position z. Note that p depends on the 



9 



difference r' ± - r\ (z'). This dependence is important concerning the effect 
studied in this paper, as it will be seen later on. 

Eq. (18) consists of two terms: one in v ± and the other in r ± - r' ± . We will 

-» 

sometimes name the first term the "current term" E ±c , while the second will 

-> 

be indicated as the gradient term E ±g . With the help of Eq. (8), Eq. (9) and 
Eq. (14) we can re-write Eq. (18) as 



z 

~ . -> . ico C , . 1 f,4 (Ke x . ,, 

E±(z, r ± ) = dz dr' ± sm(k w z ) + 

c J z-z' J \ y 
o 

xpo (r'±_ - r w cos(k w z')e x ) f(co) exp I ico 



\r ± -r' ± V 
2c(z - z') 



r± - fx 
z-z' 

icoz' ) 



ley 



(19) 



3.2 Longitudinal field 



A similar expression can be found for the longitudinal field. Since v z (z) ^ c, 
we can write the longitudinal equivalent of Eq. (5) as 



g z = -471 exp 
that is 



ICOZ 

c 



(20) 



g z = -4tt exp 



i ^ dz- 



to 



2f z c 



ico d 



^ dz 



(?± - r' 0l _ 



having used the fact that v oz (z) c. 



(21) 



It follows that the longitudinal component of the field, can be written anal- 
ogously to Eq. (18) as 



z 

E z (z, r±) = ^ dz' ^ dr' ± exp 



ico 



n -> ~t 12 



2c(z - z 



fW 



dz- 



co 



'2cy 2 z (z) 



x- 



z-z' 



ico d 

? 2 zC 



^ dz' 



Po{r\-r\ {z'))f{co) 
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(22) 

Note that the integral in Eq. (22) is performed for z' ranging from up to z 
exactly as the integral in Eq. (19), for the same reasons. 

Use of Eq. (8), Eq. (9) and Eq. (14) allow to re-write Eq. (22) as 



z 

Ez(z,rj_) = ^ dz' - — ^ dr' ± exp 



ico 



I r± - r ± \ 
2c(z - z') 



+ 



ICOZ 



ley 



r,2 



X 



ico K d 

— sinfcz )— — — 

y\c y d[x' - r' ox (z')] 



po (r f ± ~ r w cos{k w z')e x ) f(co) 



Finally, integration by parts of the term in d/{d[x' - r' ox (z')]} gives 



(23) 



z 

~ , _> , ico C , . 1 C ,-» [ 1 K . ,, 
E 2 (z, rj = rfz dr'_L — + - sm(k w z ) 

c J z-z'J [yl y 



x — x 
z-z' 



xpo (f'i. - r w cos(k w z')e x ) f(co) exp I ico 



1 r± ~ r' ± | 2 
2c(z - z') 



+ 



ICOZ 



Icy 



(24) 



Note that, in contrast with Eq. (19), we cannot clearly distinguish between 
"current" and "gradient" terms in Eq. (24): the term in sin(k w z') can be 
traced back to the gradient of the charge density, but the one in 1/yl is a 
combination between current and gradient term. 



3.3 Perturbation theory 



It is possible to analyze Eq. (19) and Eq. (24) in the framework of a perturba- 
tion theory, based on expansion in the small parameter A r /A <sc 1 according 
to the first of conditions (1). This allows simplified treatment of impedance 
and wakes. 

-> 

The first step towards this direction is a presentation of of E± and E z with 
the help of the following expansions in plane waves: 



(z - z') 



exp < ico 



r', 



2c(z - z') 



J dk ± exp [-zl ± • (r ± - r\)] exp 



ik 2 .c 

2Z7 (Z - Z) 



(25) 
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and 



r±-r' ± (. 



\?±-?± I 2 



2c(z - z') 



icr 



2nco 2 



^ dkjc^ exp J-z/c ± • {fjL - exp 
We obtain 



2Z7 (Z - Z) 



(26) 



r -» 

f ^ /M r - r r -> \Ke x . ck ± 

E±(z,r ± ) = — — dk ± dz df A — sm(k w z ) 

2ti J J J 7 w 

o L 

xpo [r'i. - r w cos(k w z')e x ) exp [-ik ± ■ (r ± - r\)] exp 



exp 



ICOZ 



2cy 



ik 2 c 



(z' - z) 



(27) 



and 



Ez(z,r_L) = J rffc_L J dz J rfrV I — - — sinfcz ) 





xpo (f'±. - r w cos{k w z')e x } exp [-ik ± ■ (f ± - r\)] exp 





icoz' 


exp 


2cf z _ 



(28) 



Performing a change of variables r' ± — > r\ - r w cos(k w z')e x and introducing 
notation 6 = k±c/co we re-write Eq. (27) and Eq. (28) as 



E±(z,r ± ) = — — — J ay J az J ar' ± I — sinfcz ) - 6 



exp 



zcuz 
2^pf 



xp (r\) exp -y • (r ± - r\ - r ro cos(k w z')e x ) 



exp 



fag 
~27~ 



(z' - z) 



(29) 



and 



E z (z, r ± ) = J 'd6 J dz' J df ± I - ^y- sin^z') 





icoz' 


exp 


_2cf z _ 



xp (r\) exp -y § • (r ± - r' ± - r w cos(k w z')e x ) 



exp 



icod 7 
2c 



(z' - z) 



•(30) 



Note that the maximal range of angles 8 X/y is limited by the last exponen- 
tial function in Eq. (29) and Eq. (30) and by the fact that z - z' > A w . It 
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follows that 6 X/ y cannot be larger than about ^A/A w . Then, the trigono- 
metric terms inside the exponential functions in both Eq. (29) and Eq. (30) 
are of magnitude cod x r w /c < ^jA r l A ^ 1- It follows that we may expand 
exp{i<jor w O x cos[k w z']/c} - 1 + iajr w 6 x cos[k w z']/c. Using exponential repre- 
sentation for all trigonometric functions we obtain: 



E±(z,r ± ) = 



cv 2 f(cv) 



ddexp 



xp (r\)exp 

icod r r 7 



2nd 2 

-8-(? ± -r> ± ) 



icod 2 z 
2c 



z 

J dz' J dr' A 



{^y [ ex P(* iyZ ') ~ expH'/^z')] - ej 



x 1 + 



2c 





~ico9 2 z'~ 




icoz' 


Jexp 


2c 


exp 


2cy\_ 



• (31) 



and 



E z (z, r ± ) = 



u 2 m 

2nd 1 



ddexp 



icod 2 z 
2c 



z 

J dz' J dr' A 



xpo (f } ±) exp -—§■ (f ± - r' ± ) | ^ - ^ [exp^z') - exp(-ik w z')] | 



X 1 + 



2c 





~ia)6 2 z'~ 




icoz' 


Jexp 


2c 


exp 


2cf z _ 



■ (32) 



Eq. (31) and Eq. (32) have been found exploiting the small parameter A r / A. 
In both Eq. (31) and Eq. (32) products of factors within {...} brackets are of 
the form exp[±ipk w z'], with p = 0,1,2..., terms for p > 1 being obtainable 
considering higher orders in ^A/A w in the previous expansion of the ex- 
ponential of trigonometric function. Note that when p = 0, the magnitude 
of 9 X/ y can be estimated from the last two exponential functions in Eq. (31) 
and Eq. (32), giving characteristic a scale d XiV ~ l/y z . For other values of p, 
instead, we have a characteristic scale Q XiV ~ ^A/A w . 

Consider first Eq. (31). Magnitudes of factors within {...} brackets are K/(2y) 
and 6 for the first bracket, 1 and cod x r w /(2c) for the second bracket. Terms 
of the form exp[±ipk w z'] with p = can have magnitudes ~ l/f z or 
[K/(2y)] ■ [co6 x r w /(2c)] ~ K 2 A r y z /(y 2 A), this last kind being negligible. When 
p = 1 terms have magnitudes K/(2y) or 9cod x r w /(2c) ~ K/y, and both kinds 
have to be kept. Similarly, it can be shown that all other values of p give 
negligible terms. 



Consider now Eq. (32). Magnitudes of factors within {...} brackets are l/y 2 
and Kd x /(2y) for the first bracket, 1 and cod x r w /(2c) for the second bracket. 
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Terms of the form exp[±ipk w z'] with p = can have magnitudes l/yl or 
[Kd x /(2y)] ■ [cv6 x r w /(2c)] ~ K 2 A r /(y 2 A), this last kin d bein g negligible. When 
p = 1 terms have magnitudes K0 x /2y ~ K/y ■ -\/A/A w or cvd x r w /(2cy 2 ) ~ 
VA r /A • l/(y z y), and this last kind can be neglected. Similarly, it can be 
shown that all other values of p give negligible terms. 

Altogether, we obtain the following expressions for transverse and longitu- 
dinal field: 



£±(z,r±) = 



2nd 1 



z 

jdd j dz' jdr\ 

o 

xpo (f'i) exp -™§ • (r ± - r\)] exp 



\ia)6 2 (z'-z) 




icoz' 


[ 2c \ 


exp 


2cy\_ 



xl-d + 



Ke x _ icod x r w d 
2iy 2c 



exp(ik w z') - 



Ke x icod x r w 6 
2iy 2c 



exp(-ik w z') 



(33) 



and 



E z (z, r±) = 
x exp 



u 2 m 

2nd 2 

icv6 2 (z' - z) 
2c 



z 



exp 



icoz' 111 K6 X 



2cy 2 z \ [y 2 2yi 



[exp(ik w z') - exp(-ik w z')] | . 

(34) 



Eq. (33) and Eq. (34) are the first order result of our perturbation theory, 
where the small parameter A r / A has been exploited through the expansion 
of exponential functions in Eq. (29) and Eq. (30), and non-negligible terms 
are kept. 

We now go back to the space-frequency domain performing the integral in 
dd with the help of Eq. (25), Eq. (26) and using also 



J deexp[-^e-(f ± -r\) 

2cn [c(z - z') + ico(x - x') 2 ] 
= (z' - z) 3 

Finally, we obtain: 





ico9 2 (z' -z) 


exp 


2c 



exp 



ico \r ± - r' ± \ 
2c(z - z') 



el 



(35) 
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E±(z,r ± ) = - 



icofico) 



c 

z 



J dr'±p (r\) 



exp 



ICOZ 



X- 



exp[+ft „, 2] j _ exp 

o 

-exp[-ik w z] J ~~^ ex P 

o 

- exp[+z'fc l( ,z] I _ - exp 



2cy 

2 



r± - ?" ± r 

2c(z - z') 



Ke x 

+^exp[ik w (z' -z)] 
Ke 

-^exp[z'fc ro (z-z')] 



x - 



IjWx I 2 
2c(z - z') 



2(z - z') 2c(z - z') 2 



exp[z'fc n ,(z' - z)] 



expH'fc^z] J" ^^exp 



ico- 



r ± - r' ± Y 
2c(z - z') 



+ 



J z-z' 



X 



exp 



2(z - z') 



2c(z - z') 



>\2 



exp\ik w (z - z')] 



r i - r 



zcu- 



2c(z - z') 



exp 



z<i>(z' - z) 
2cy z 2 



z-z' 



(36) 



for the transverse field, and 



E z (z, r ± ) = 



icof(co) 



c 

z 



J dr\p (fj) 



exp 



zcuz 



2cf z 



r dz' 

exp[+z'fc z „z] I j—^ exp 



ico- 



.\ ?±-?± P 
2c(z - z') 



+ - K % X , exp\ik w (z' - z)] 



2z'y z-z' 



+ 



+ 



expH'/c^z] J ^jp-exp 

o 

J z~^I' 



ico- 



r ± - r' ± \ z 
2c(z - z') 



K x-x' 



2z'y z-z 



- exp[-ik w (z' - z)] 



r , - 1 



exp \tco- 



2c(z - z') 



exp 



za>(z' - z) 
2cy 2 



y z 2 . 



(37) 



for the longitudinal field. Here we neglected factors exp[ico(z' - z)/ (2cy 2 )] in 
integral terms in dz' including exp[±k w z], because coA w /(2cyl) <sc 1. 

Note that there exists a mathematical shortcut to obtain Eq. (36) and Eq. (37) 
from Eq. (19) and Eq. (24). In fact, if we perform a change of variables r' ± — > 
r'± - r w cos(k w z')e x , we formally expand the Green's function exponential 
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exp{z'<y[| r ± - r' ± - r w cos(k w z')e x \ 2 /2c(z - z')]} to the first order in r w and keep 
non-negligible first-harmonic terms in exp[±ik w z'], we obtain Eq. (36) and 
Eq. (37). We will regard it as a mnemonic rule, that will be useful later on. 

Fields are calculated under conditions (1). In the limit for z » y\K, as we 
will see in the next Section 3.4, integrals in dz' depend on z only through 
phase factors, i.e. a steady state solution is reached. 

Analysis of Eq. (36) and Eq. (37) presents an interesting picture of the fields 
generated by the electron beam. Eq. (36) and Eq. (37) consist of the sum of 
integrals in dz' . Some include exponential factors exp[±k w z], other not. 

Terms not including exp[±k w z] (the last integrals in dz' in both Eq. (36) and 

■± 

Eq. (37)) oscillate, as a function of z, on a scale Ky\. The field E is given 

■± 

by Eexp[z'cuz/c]. It follows that the phase velocity of terms not including 
exp[±k w z] is the same as that of the electron beam harmonic p. We can 
interpret this fact by saying that this part of the field is entangled with the 
electron beam. It is natural to identify these terms as space-charge terms. 
The formation length of the space-charge field is determined by the factor 
exp[z<y(z' - z)/(2cfl)] under integral sign, and amounts to 2.Xy\. Similarly, 
the diffraction size of the space-charge field is given by y z K. 

Terms including exp[±ik w z] are indicative of fields E performing a cycle of 
oscillation on the scale of an undulator period with respect to the electron- 
beam harmonic p. Phase velocity of terms including exp[+ik w z] is slower 
than that of the beam harmonic. These field terms have a phase velocity 
slower than the speed of light. Phase velocity of terms including exp[-ik w z] 
is faster than that of the beam harmonic. These field terms have a phase 
velocity faster than the speed of light. We can interpret these facts by saying 
that these parts of the field are not entangled with the electron beam. It 
is natural to identify these terms as radiation terms. The formation length 
of radiation field terms is determined by the factor exp[±ik w (z' - z)] under 
integral sign, and amounts to K w . Similarly, the diffraction size of the space- 
charge field is given by yjKK w . 

It is interesting to trace each term in Eq. (36) and Eq. (37) back to the source 
terms that originated them, distinguishing between gradient and current 
terms. Considering Eq. (36) it can be seen that the first and the second 
integral are (radiative) current density terms. The third and the fourth term 
are (radiative) gradient terms, while the last term is a (space charge) gradient 
term. However, one can see from Eq. (37) that the third integral is a (space- 
charge) term originated from a mixture of gradient and current sources. 
Thus, although the first and the second integral are (radiative) gradient 
terms, it does not make sense to separately talk about gradient and current 
term for the longitudinal component of the field. 
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An interesting picture emerges, where radiation field and space-charge field 
are treated on equal foot, through paraxial Maxwell's equation. On the one 
hand, as we have seen, these fields have different formation lengths, and 
different diffraction sizes. On the other hand, our theory allows for generic 
transverse sizes of the electron beam Oj_, that makes it possible to compare 
± with both diffraction sizes, thus obtaining different regimes. As we will 
see, when o ± » ^KK Wl impedance and wakes are essentially dominated by 
the longitudinal space-charge term. It is important to remark, for future use 
in the next Sections, that y z enters the expression of the space-charge field, 
and not y. 



3.4 Explicit expressions for the field 



We can consider Eq. (36) and Eq. (37) as starting point for our investigations, 
and calculate explicit expressions for the field to be used later on in the 
calculation of the impedance. First, we make a change in the integration 
variable from z' to £ = z - z' . In the limit for z — > oo, corresponding to the 
second of conditions (1), i.e. z » y\K, we can write 



E±(z,?±) = - 



icofioS) 



exp 



10)Z 
2 

z . 



ley 



Ke x . „ . Cdl 
+ exp[+ik w z] J — exp 

o 

oo 

Ke x Cdt 
- expL-z^z] J — exp 



. r ± - r\ | 2 
+uo z—: ik w t, 



+ico - 



2cl 



r, -f 



+ exp[+ik w z] 
. 2icr w 



+ 



2cl 
d 



0) 



2co\A-r\\ d[\? ± -r> ± \] 

rP- 

(x - x )(r ± - r\) • — - 7 

d [| ? ± - r' ± | 2 ] 



+ exp[-ik w z] 

2icr w 



x 



+ 



dl 




-* ~i 12 
r ± - r'j_ f 


— exp 


+io)- 




-> 

icr w e x 




d 



ikw£ 



+ 



0) 



2o)\? ± -r ± \ d[\r ± -r' ± \] 

d 2 

(x - x')(r ± - r\) ■ — - 7 

d[\n-r' ± | 2 ] 
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X 



oo 

If 



exp 



+ I(0~~~z — ~~ h 



2c£ 



+ 



ic r ± - r' ± 



|r ± -r\| i[|rl-r\|] 

oo 

x j T exp 



j± - r'± I 2 



2c£ 2cy 



v?2 

z . 



(38) 



for the transverse field and 



E z (z, r ± ) = 



icof(co) 



J df ± p (r\) 



x| - exp[zX w z] 

oo 

J T exp 

o 

- exp[-^„z] 

oo 

x j -ex, 

o 

oo 

y 2 J ^ 



exp 
x — x' 



ICOZ 



2coy\f ± - r -! ± \ d[\? ± -r> ± \] 
. \r ± -f ± \ 2 



ICO 



2c£ 
cK x — x' 



+ICO- 



2cl 



+100 



r ± - r' ± | 2 iooE, 



2cl Icy 



for the longitudinal field. 

We now use the fact that, for any real number a > 0: 



(39) 



J °° d£ exp [i + a/O] It = 2K (l yfa) 

J °° &l exp [i (+£ + a/O] /I = 2K (-2i yfa) = n [ij (2 y/a) - Y (2 y/a) 



(40) 



where K n is the n-th order modified Bessel function of the second kind, Y n is 
n-th order Bessel function of the second kind and /„ is the n-th order Bessel 
function of the first kind. Using Eq. (40) and the fact that k w » co/(2cfl), and 
remembering that 
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d 2 

d[\r ± -r\ |f 

1 ' d 2 



4 | ? ± - r\ \ 2 d [| ? ± _ r ->, ± |] 2 4 | ? ± - r\ f d [| ^ - r \ |] 
we can write Eq. (38) and Eq. (39) as 



E±(z,r ± ) = — 



)exp 



zcuz 



Icy 



I *y \ Vaa 



- expH'/c^z] 1 - V2z . 

iy \ Vaa 



w 

-> 



w 



+ exp[+ik w z] 



icr 7t 



exp[-ik w z] 



CO 



ten, 



y/2'e 



-Ki V2- 



r, -r'. 



VAA^ \r±—r' 

+ 2 ( x-x-)y ± - r \) K /^ 



Vaa^ 

r , - r' 1 



AAa, I r ± - r' ± | 2 
VAA^ I r ± - r\ I 



Ki -V2z 



VAA ro 
\r 1 - r' 1 



Vaa; 



+ 



2(x-s')(r ± -r' ± ) ^ /_ ^. |r ± -r' ± | \ 
AA ro I r ± - r*_L | 2 I VAA^ / 



z'c r ± - rV 
°> I r ± - r\ I 



and 



(41) 



(42) 



E z (z, r ± ) = - 



icof(co) 



J df^po (r\) 



exp 



zcuz 



2cyl 



V2 

x<( h — exp[+ik w z] 



+ 



VAAa, 

V2z 

VAA Z „ 



exp[-ik m z] 



cK x-x' 

W Y I r ± - r' ± 
cK x-x' 
<vy\r ± - r\ 







r ± - r' ± \\ 


VAA ro / 






r ± -r' ± \\ 


VAAa, / 



Yz \ 7zA 



(43) 



Note that, similarly as in Eq. (36) Eq. (37), it is possible to recognize in Eq. 
(42) and Eq. (43) radiative and space-charge terms, as well as gradient and 
current terms (for the transverse field components). 
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3.5 Cross-check with Gauss law 

It is possible to cross-check our expressions for the field with the help of 
Gauss law: 



This cross-check will constitute a general cross-check of the correctness of 
our calculation and allow a better understanding of the interplay of different 
field contributions in the complicated machinery of Maxwell's equation. 

As said before, in the present study we work in the steady state, when the 
bunch has travelled inside the undulator for more than 1y\X. In this case, 
an explicit expression for transverse and longitudinal fields are given in Eq. 
(42) and Eq. (43). 

We will demonstrate that the field E = E ia< x + E sc obeys Eq. (44) by separately 
showing that radiation field E ra d and space-charge field E sc verify: 



Relations (45) can be interpreted saying that the radiation field is not entan- 
gled with sources, while the space-charge field is. Hence the different right 
hand sides. 

Let us begin with the space-charge field. First we can write: 



because, in the steady state, E z sc depends on z only through exp[z'cuz/ (2cyl]. 
Thus, in order to verify the first of Eq. (45) we should prove that 



V • E = 4np . 



(44) 



V • E sc = 47ip 

V • E rad = . 



(45) 




(46) 



V ± • E_l sc = 47ip £ zsc . 



(47) 



From Eq. (42) we have 
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F 



2f(co) 



exp 



ICOZ 



2cy 



J r V ' | ? ± - r\ | \ 7zA 



(48) 



In order to calculate the left hand side of Eq. (47) we can use the divergence 
theorem in two dimensions, to find: 



V, • 



CO 



Cfz 
V ) 



CO 



R 



cy z 

V ) 



+ 



2ncy 2 



CO 



6(Rx) , 



(49) 



where 5 indicates the Dirac-delta function and derivation is understood in 

weak sense, we set R ± = r± - r'±. Remembering E± = E± exp[-icoz/c] we 
obtain 



—5-; dr\p (r' ± )K 

c 2 yi J 



(\-> ~i i\ 



Ay 2 



+ 471 p (f ± ) 



~xj(co) exp 

Now, from Eq. (43) we have: 



(50) 



F = 



2z'cu 



exp 



zcuz 



2cy; 



K, 1 ! ' r± - r ± ] f(co) exp 



Ay 2 



zcuz 



L V z J 



(51) 



Substitution in the right hand side of Eq. (47) yields Eq. (50), thus verifying 
Hq- (47). 

Let us now consider the radiative fields, and show that the second of Eq. 
(45) is also verified. Presentations in Eq. (42) and Eq. (43) include now many 
terms, and it is convenient to start with alternative presentations for the 
transverse and longitudinal field, namely Eq. (33) and Eq. (34). In the limit 
for z » y\K The radiative part of the field is given by: 



E±rad(z, r ± ) = exp 



ICOZ 



u 2 m 

2nc 2 



00 

j ' d§ J dz' J df ± p (r\) 



x exp 



ico -+ /_, -> \ 



exp 



icod\z' - z) 
2c 
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X <^ 



Ke x . 

— sm[k w z J - 



ico6 r r m 6 



cos[k w z'] 



(52) 



and 



J z rad 



x exp 



(z, r ± ) = - 



co 2 f(co)K 
2nc 2 y 



exp 



GO 

~c~ f d ® f dZ ' f d ^ ±p ° ( r1± ) 



-6 ■ (r ± - r> ± ) 



exp 



!&)0 2 (Z' - Z) 

2c 



6 X sm[k w z'\ . 



(53) 



Here we neglected factors exp[/cuz' / (2cyl)] because coA w / (2cy 2 ) <sc 1. We may 
now calculate directly V ± • E ±rad and d z E z rad . We obtain 



V ± • E ±rad (z, r ± ) = - 



ico Kco 2 f(co) 
c 2nc 2 y 



exp 



ICOZ 



GO 

Jdd J dz' J dr' ± p (f' ± ) 



x exp 



—-e-(r ± -r> ± ) 



-26 X A W cos[k w z']-^ exp 



d x sin[^„z'] exp 

ico0 2 {z' -z) 
2c 



~ icod 2 (z' -z) 
2c 



(54) 



and 



ico Kco 2 f(co) 



c 2nc 2 y 



icoz 



d z E z rad(z, r±) = — exp 

~ico0 2 {z' -z) 



CO 

Jdd J dz' J dr' ± p (r\) 



x exp 



ICO ^ /_, -> \ 

-e ■ (f ± - f ± ) 



exp 



2c 



{d x sin[k w z']} , 



(55) 



In Eq. (55) we neglected an extra-term -6 2 6 x /2 sm[k w z'] in parenthesis {...}, 
because A/A w «: 1. We now integrate by parts the term in cos^z'] in Eq. 
(54) to obtain 



V ± • Exrad(Z, T±) = - 



ico Kco 2 f(co) 
c 2nc 2 y 



exp 



ICOZ 



oo 

Jdd J dz' J dr' ± p (f' ± ) 



x exp 



ICO £ /_> -> \ 



exp 



icod 2 (z' - z) 
2c 



exp 



o 

icoz' 
2cf z 



{-^sin^z']} . 

(56) 



Obviously V ± • E ±r ad + d z E z rad = 0, and also the second of Eq. (45) is satisfied. 
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4 Wake function and impedance 



The knowledge of the electric field at the position of any test particle inside 
the beam, that has been derived in Section (3), allows to derive wake function 
and impedance related with the system under study Let us briefly review 
here concepts of longitudinal impedance and wake function, that will be 
used later on. The longitudinal impedance of a system, Z (co), can be given 
as the Fourier transform of the wake function G (As): 



CO 



d(As) 



G (As) exp 



ico- 



As 
i — 



G„ 



d?- E°(As,r'(t),t)\ t = z , /iM . 



(57) 



Here the integral in the expression for G is a line integral calculated along 
the trajectory of a test particle. In fact, E°{As,r'{t),t) indicates the time- 
domain electric field generated by a source particle acting on the test particle 
at longitudinal distance As from the source. In the calculation of E° (As, r'(t), t) 
we assume that effects from the vacuum chamber are negligible, i.e. a » y z A, 
that is the third of conditions (1). E°(As,r'(t),t) is integrated along the test 
particle trajectory, and divided by the electron charge (— e), so that e 2 G (As) 
is the energy (gained, or lost) by the test particle due to the action of the 
source. In agreement with |20|1 we take the test particle behind the source 
for positive values of As. 

According to the given definition of wake function, one should integrate 
the field over the entire trajectory. However, there is no principle difficulty 
in considering only part of the trajectory, let us say, up to longitudinal 
position z. Mathematically, this means that the line integral for G should 
be performed up to the trajectory point of the test electron corresponding 
to longitudinal position z. In this way, G = G(As,z). 

Note that Eq. (57) is automatically dependent on a particular source electron, 
and a particular test electron. In order to formulate this statement in a math- 
ematical way, we may introduce test and source particle initial transverse 
position r ±T and r ±s and write G = G (As,z,rj_s,f ±T ), where we neglect 
differences in energy between the two particles. Following our previous 



work 121(], we will slightly modify the concepts of wake and impedance by 
substituting test and source particles with disks of total charge (-e), longitu- 
dinally separated by a distance As. This amounts to an integration over the 
transverse particle distribution in dr ± T and dr ± s, that makes our definitions 
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independent of r^r and r ±S - We thus obtain 



G(As, z) = c 2 j dr\ j dr' ± p (? ± )p (/' ± )G (As, z, r\, r" ± ) 



(58) 



In Eq. (58) we used the fact that p is independent of the longitudinal posi- 
tion. With the redefinition in Eq. (58) we can further consider the impedance 
Z(co,z) proceeding as in (57) for the definition of Z (co), but Fourier substi- 
tuting G (As) with G(As,z). Note that by definition of Z(co,z) we have 



Z(co, 



,z) = -J^ — f T -EdV = -J- — f cfe' \ dr'X -E, 



(59) 



where |/(<y)| 2 accounts for the fact that test and source disks have total 
-±* -± 

charge (-e), while j ■ E oc |/(cu)| 2 , f(co) being the already defined Fourier 
transform of the longitudinal bunch-profile. Here the volume V is a cylinder 
of base A including the undulator up to position z' = z. The integration in 
z' is performed from to z, because we will be interested in impedance 
and wakes generated inside the undulator, and we will assume that the 
undulator begins at position z = 0. 



5 Impedance calculations 

According to Eq. (59), the expressions for the longitudinal impedance asso- 
ciated with E ± and E z , Z ± and Z z are given by 



z 

Z ± = ^— fdz' 

l/MI 2 J 



exp 



is (z') 



x J dr\ f*(co)p* [r' ± - r w cos(k w z)]v 0± (z') ■ E ± (z', co, r' ± ) 

i r r i So (z'Y 

= — I dz exp 

|/(o;)| 2 J F [ v 



(60) 
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When we calculated the field, we saw that the dependence of the source 
transverse charge density p on the electron motion r' ±0 (z) had to be ac- 
counted for. We will show that, in the calculation of the impedance, the de- 
pendence of the test transverse charge density p on f ' ± - r w cos(k w z)e x must 
also be accounted for. With the help of a change of variables we rewrite Eq. 
(60) as a sum Z = Z ± + Z z : 



z 

Z = ^— fdz'exp [ dr'J^ooWAr'J 



{Kc 
- - [exp( ft „z<) - exp(- ft „z')] B#, a . r\ + ,.e»fcz)?,) 

+ cE z (z', oo, r\ + r w cos(k w z)e x ) \ . 



(61) 



Calculations can be drastically simplified, because transverse radiative gra- 
dient terms in Z ± cancel with longitudinal radiative terms in Z z . It is easier 
to show this facts with the help of Eq. (36) and Eq. (37), rather than using 
explicit expressions Eq. (42) and Eq. (43). 

First, with the help of Eq. (36) and Eq. (37), we write down the part of the 
impedance from the radiative transverse field, Z ±r and from the radiative 
longitudinal field, Z z r . In principle, in order to dispose of the oscillating 

terms in E T (z',co, r\ + r w cos(k w z)e x ), we may use the same mathematical 
shortcut that can be exploited to obtain Eq. (36) and Eq. (37) from Eq. (19) and 
Eq. (24). In fact, we may formally expand the Green's function exponential 
exp{zcu[| r' ± - r" ± + r w cos(k w z')e x \ 2 /2c(z' - z")]} to the first order in r w and 

keep non-negligible first-harmonic terms in exp[±ik w z']. However, E ±r is 
multiplied by v OJ _ in the expression for the impedance Z ± . Since v OJ _ oscillates 
with period A w , we can neglect the oscillatory contributions in cos(k w z)e x in 

the expression of Ej_ r , because they would give oscillatory contributions that 
average to zero after integration in dz' . We therefore obtain: 



z 

Z ±r = |^ J dr' ± j dr" ± p* (r' ± )p (f' ± ) J dz' 



exp \ ioo- 



f i - r 



± i ± 



2c(z' - z") 



\ Ke 

+ _Zexp[^(z"-z')] 
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+ 
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icor w (x' - x") 1 
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exp[zfc a ,(z' - z")] 



(62) 



for the transverse field. On the contrary, the longitudinal velocity v z is a sum 
of a constant term, whose magnitude is about c and a negligible oscillates 
with period 2A W . As a result, the oscillatory contributions in cos(fca,z)^; in 
the expression of E z r must be kept, and an expansion of the exponential in 
the Green's function must be performed, leading to 
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exp zcu— — - 
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exp[zla,(z" - z')] 
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exp 
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Z6jr It; (x' - x"f _ r m 
2c(z' - z") 2 2(z' - z") 



exp[zUz' - z")] 



(63) 



Partial cancellation can be exploited between the last two terms of Eq. (62) 
and Eq. (63). Note that the longitudinal radiative impedance is completely 
cancelled, and one obtains 
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exp 



exp 
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zcu- 
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+—exp[ik w (z" -z')] 
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(64) 



Performing the integrals in dz" with the help of Eq. (40) or, equivalently, 
calculating Z r with the help of the transverse radiative current terms in Eq. 
(42) we finally obtain 



K 2 coz 
2y 2 

V2 | r\ - r" ± | 




J J" df\p* {f Jpo (r" ± ) 



- V2z | r' ± - r" 

VAA ro 



(65) 



Let us now consider the space-charge part of the impedance. One finds that 
the space-charge term in Z ±/ i.e. Zj_ sc , averages to zero, as can directly be 
seen by inspecting the last term in Ki [|r ± - r' ± |/(y 2 A)] of Eq. (42). In fact, such 
term is independent of z. Now, according to Eq. (61), in order to obtain the 
correspondent impedance contribution, this term must be multiplied by v± 
(i.e. by exp[±ik w z']) and integrated in dz' for a saturation length. It follows 
that, during the integration process in dz' , one integrates a fast varying 
function of z on the scale of K lv . As a result, we obtain a negligible effective 
impedance contribution over many undulator periods, and we can neglect 
the transverse space-charge contribution in the calculation of wake and 
impedance. The total space-charge part of the impedance coincides with 
the longitudinal space-charge impedance. It can be shown that oscillatory 

contributions in cos(k w z')e x in the expression of E 2SC are of higher order in 
co I co r and they can thus be neglected. As a result we obtain 



__. 2 + K 2 

^sc — sc — 10JZ 



j df ± J dr'\p:{r\)p o (f\)K ^— 



_L ' _L 



(66) 

We thus reach the conclusion that only longitudinal space-charge terms and 
transverse radiative terms enter the expression for the impedance that can 
now be calculated in integral from for any transverse beam-distribution 
under conditions (1). 
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Finally, we obtain the total impedance Z = Z r + Z sc . The real part Z R is given 
by 



Z R = - 



K 2 



71 CUZ 



4y 2 



J df ± j dr" ± p* (F> \)p (r" ± )/ 



VAA l( , 



The imaginary part Zj, instead, amounts to 



Zz = - 



2y 2 
-Ko 



J~ dr\ J* dr*' J _p' (ri A _)p (r" ± )j^Y 
V2 |r\ - r" ± \ ) 4 + 2K 2 f Ir 1 , - r" 1 N 



V2|r\ 



VAA- 



VAAj, 



K 2 



y z A 



(67) 



(68) 



having used the fact that K (-ix) = (n/2)[ij (x) - Y (x)]. 



5.1 Asymptotic case for cr 2 « AA Z , 



Before proceeding with the analysis of the wake, it is interesting to derive 
asymptotic limits of Eq. (67) and Eq. (68) in the case for cr 2 «: AA W . Bessel 
functions in Eq. (67) and Eq. (68) can be expanded for small argument values. 
In particular, using } (x) 1 for x <sc 1, and recalling that p is normalized 
to 1/c, the real part of the impedance becomes 



Z R = - 



K 2 conz 
4y 2 



J dr^ J dr" ± p* (r' ± )p (f' ± ) = 



K 2 nz 
'4c Ay 2 ' 



(69) 



independently of the choice of p . Subsequently, we use K (x) -y E -ln(x/2) 
and Y - 2/7i[y E + ln(x/2)], y E - 0.577216 being the Euler Gamma constant 
in the imaginary part of the impedance, Eq. (68). We obtain 
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where 



_ 2zy E 2coz 

free — Z 9 "I n~ 

cAy z y z 



/ 1 -> -> 

T 1 - f' 



2Ay 



(71) 



Zj free is the only model-part of the impedance. In particular, assuming a 
Gaussian transverse profile: 



po(r±) = 



2no 2 c 



exp 



2a 2 



(72) 



we obtain 



2zy E z 
= — - — 1 

cAy 2 In^cKy 2 



x 




2zy E | 2z ln /^x\ 
cAy 2 cAy 2 \yA/ 



(73) 

Also, Z; f re e, is logarithmically divergent on cr ± . This is, in fact, the free-space 
impedance. The renormalized impedanceEl i.e. the difference of Eq. (70) 
with the free-space impedance is independent of a ± and constitute a result 
valid for any value of K. 

Here we underline the fact that in the limit for o ± — > the difference 
between the impedance of an electron beam moving through a magnetic 
system and the free-space impedance is independent of the electron beam 
model. It is finite, and can thus be applied in a one-dimensional approx- 
imation whereby the electron bunch is modelled by a line density. This 
one-dimensional approach was first proposed in the time-domain to study 
Coherent Synchrotron Radiation (CSR) in [221 ] and is currently used in CSR 
codes. Also in this case, the renormalized wake is obtained by subtracting 
the free-space wake (as is natural, because the wake is the Fourier transform 

5 Note that, without slow-wave radiative contributions to the field (proportional to 
exp[+ik w z]), it would be impossible to recover Eq. (70) and Eq. (80). In other words, 
the renormalization process would fail. This underlines the fact that, although 
slow-wave radiative contributions have no realization in the far zone, they are of 
fundamental importance in the calculation of the impedance. 
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of the impedance, i.e. its time-domain counterpart). The renormalizatiorfll 
procedure used here and introduced in | |22|l has to be seen as a mathematical 
algorithm to deal with calculation of self-forces of a moving charged line. 
Such calculation is problematic due to incompleteness of electromagnetic 
theory, yielding to divergence. Such divergence is cancelled by subtract- 
ing the longitudinal force that would be present in a straight-line motion 
from the force calculated on a curved trajectory The finite difference can be 
entirely ascribed to curvature. 



5.2 Existing studies of the asymptote o\ «: AA W 



Analytical and numerical studies can be found in literature, treating longi- 
tudinal impedance and wake from an undulator setup in the case of a line 



density distribution of electrons (see e.g. |24], [25|] and [26]) 



Reference |24J1 deals with the one-dimensional renormalized wake of an 



electron beam with a Gaussian longitudinal profile 



„v HON 

f(s) = exp 

V27IC7 Z 



(74) 



N being the number of electrons in the beam. In particular, in reference 



||24j], the following expression for the energy gained or lost by a particle 
at position z down the beamline and position s within the buncrQ was 
obtained: 



AS = J_ G(p, K, x) , (75) 



V2~7zcr z y 2 

where x = —s/o z , p » 1 is the bunch length parameter 



v = TTWJ2 (76) 



and G is given by 



6 It may be worth to note here, that in the renormalization procedure used in [22 
only retarded fields are used whereas in previous works (see e.g. [23]) devoted to 
renormalization in classical electrodynamics a radiation field is used, that is half 
the difference of retarded and advanced fields. 

7 It should be noted that the definition of s in this paper differs for a sign with 
respect to that in [24]. 
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G(p,K,x) = -exp 



XT 

T 



[ln(p) + g(K)] + F(x) , 



(77) 



with 



F(x) = ^[y E + 31n(2)- 


- 2]x exp 


' X 2 

2 




-VIH 


X 

L V2 J 


- xexp 


X 2 ' 

2 



X 



exp 



^2 



1 + erf 



x' 



LV2 



(78) 



Moreover, in the limit for K 2 «: 1, g(K) — > 0, while in the limit for K 2 » 1, 
g(K) — > 1. For arbitrary values of K, g(K) was presented in [2J] as a plot, 
using numerical integration techniques. Now g(K) can be expressed fully 
analytically: 



g(K) = 1 



In [1 + K 2 /2] 
W/2 



(79) 



Eq. (75) was already been independently cross-checked, with the help of the 
code TraFiC 4 , in [2jft. In this paper we underline the correctness of Eq. (75). 



Reference |26|1 deals with the renormalized impedance in the case of a line 
density distribution. Such impedance is presented in Eq. (26) of |26|l in the 
asymptotic case for K » 1. When o ± «: ylKK w and K 2 » 1 only the first 
term of Eq. (70) survives, and the total renormalized impedance Z ren reads: 



„ _ K 2 nz . K 2 z I A 

Zj — IZjI f ree — ; — ; — - — I— — ; — - in 



4cAy 2 2cKy 2 



K r ) 



(80) 



Eq. (80) is in agreement with reference |26|], where the impedance per unit 
length is givero. 



5.3 Energy conservation law for a]_ «: AA n 



In general, the real part of the impedance can always be cross-checked with 
the energy conservation law, that requires: 



8 An extra factor — 1/c in our expression is the result of different definition of 
impedance. 
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^ = -^fipf Z R (fi>) , (81) 

where the energy spectrum of the radiation, dW/dco, is defined as the integral 
over all angles of the total energy emitted per unit frequency per unit solid 
angle dQ. = ddddcp : 



dW 

dco 



2n 



2n 



de e 



(82) 



It is easy to verify Eq. (81) in the case a]_ AA W . In this case, the electron 
beam transverse size is much smaller than the radiation diffraction size, and 
a filament-beam model can be used. 

Our theory has been developed for any value of the undulator parameter 
K, and in the long-wavelength asymptote, i.e. for A r /A «: 1. In this case it is 
possible to give a simple mathematical description of the radiation energy 
spectrum. Based on this expression we will then verify Eq. (81). 

In order to calculate the energy spectrum according to Eq. (82) we must first 
-> 

calculate E ± in the far zone. We can specify "how near" co is to the resonant 
frequency co r = 2k w cy\ by introducing a detuning parameter C, defined as 
C = co I \2y 2 c) —k w = (Aco /co r )k w , where co = co r +Aco. Then, the field generated 
by a filament-beam is well-known and is given, in paraxial approximation 
(see e.g. Eq. (13) of Q), by: 



E ± = exp 



.cod 2 z 



ico\f(co)\ 



Lw/2 



x exp 



i\C + 



2c 


C 2 Z 






-1 


W2 


cod 2 ' 




CO 


2c , 




k w c 



cos(k w z') 



8y 2 k w c 



sm(2k lu z') 



(83) 



As first proposed in |27I1 one may use the Anger-Jacobi expansion exp [ia sin(i^)] 
Tj^L-co Jp( a ) exp [ipip], where } p (-) indicates the Bessel function of the first kind 
of order p, to write the integral in Eq. (83) in a different way: 



exp 



.co 
i- 



e 2 z 



2c 



ico\f(co)\ 



C 2 Z 



Y^j Jm(u)Jn(v) exp 



m,n=-oo 



inn 



Lw/2 

x ^ dz' exp 



i\C + —— lz 



2c / J [2iy 



K 



— [exp (2ik w z') - 1] e x + exp (ik w z') 
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x exp [i(n + 2m)k w z'] , 



(84) 



where u = -K 2 co/(8y 2 k w c) and v = -KO x co/(yk w c). There is no simple result 
valid at the same time for arbitrary K values and arbitrary detuning. How- 
ever, there are asymptotes for A r /A «: 1 and arbitrary K, or K 2 « 1 and 
arbitrary detuning. We will first assume K 2 « 1 and arbitrary values for 
A/A,0- Consider the case \Aco\/co r » 1/N W , N w being the number of undu- 
lator periods, C < and K 2 «: 1. Because of these assumptions, both u «: 1 
and u« 1 (here 2 ~ 2c\C\/co). This means that asymptotic expansions of 
Bessel functions can be used. Non-negligible contributions are for n = m = 
or for n = -1 and m = 0. It follows that 



Ei =■ 



co\f\co)\K 
2c 2 zy 

L w /2 

: J dz' < 



exp 



.cod 2 z 



e x -6 



6 x co 



exp 

-Lw/2 

2c 2 zy 
L 



xsmc 



exp 
4 2c 



.co0 2 z 
i 

2c 



kjnC 



e x e y o)' 

k m c 



(85) 



Note that since both u <s 1 and p « 1, the Anger-Jacobi expansion is not 
really necessary here, and we might have derived Eq. (85), based on K 2 «: 1 
and \Aco\/co r » 1/N W , directly from Eq. (83) by directly expanding in Taylor 
series the exponential function of trigonometric arguments. 

The total energy emitted per unit frequency per unit solid angle K 2 <s 1 and 

\Aco\/co r » 1/N W is 



o; 2 |/(o;)| 2 K 2 L 



dcodD. 16n 2 c 3 y 2 



k w c 



+ 



k w c 



smc 



(86) 



in agreement with [28]. Substituting Eq. (86) in Eq. (82) and using the fact 
that N w » 1 and sinc 2 [x/fl]/(7Zfl) — > 5(x) for a — > we obtain 



9 Note that in the long wavelength asymptote, A » A r , i.e. cu/cy r «c 1, we always 
have \A(o\/(o r » 1/N W , but not viceversa. Thus, here we are considering \A(o\/co r » 
1/N WJ but arbitrary value for A/A r . 
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dW co\f(co)\ 2 K 2 L 



dco 8c 2 y 2 



1 + 



CO 



x ck w y 



(87) 



Also Eq. (87) is valid for K «c 1 and \Aco\/co r » l/N w , and is in agreement 



with 1270 (where the energy spectrum was first calculated) and, | |28|| 1 and 



Q24f [. In the limit for A » A,- we write 

dW co\f(co)\ 2 K 2 L w 



dco 4c 2 y 2 



(88) 



Eq. (88) is at the left hand side of Eq. (81) and has been calculated for K «: 1 
and A » A,-. 

The right hand side can be written using the real part of Eq. (69), that is valid 
for arbitrary K values and A » A r , and calculating the impedance along an 
undulator of length L w . We obtain: 



1 . |2 dZ R co\f(co)\K 2 L w 
~ n Lw Hz- = 4c 2 y 2 ' (89) 

thus verifying Eq. (81) and energy conservation. Agreement between the 
asymptote of Eq. (87) for A » K r and Eq. (69) is due to the fact that Eq. (69) 
is valid for arbitrary values of K. However, Eq. (88) can also be calculated 
from Eq. (83) under the only assumption A » A,-, i.e. Eq. (88) is not only 
valid for K 2 <s 1, but for any value of K. More generally, we can say that in 
the long wavelength asymptote (A » A r ) it is sufficient to account for the 
first harmonic only, independently of the value of K. 



5.4 Asymptotic case for o\ » AA Z , 



As we already discussed, radiation field and space-charge field exhibit dif- 
ferent formation lengths and different transverse scales, namely VAA IW and 
Ay 2 . The same transverse scales are also present in Eq. (67) and Eq. (68). 

The first two terms in Y and K in Eq. (68), as well as the entire real part of the 
impedance, are linked to the presence of transverse current density and to 
radiation field. The last term in Eq. (68) instead, is due to the presence of lon- 
gitudinal space-charge field, a combination of current and gradient terms. 
The corresponding Bessel functions yield different characteristic transverse 
scales. Bessel functions related with the radiation field are linked with a 
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A typing error is present in Eq. (2.11) of [28]. 
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transverse size yjAA w ~ ^AA r y\. Those related with the longitudinal space- 
charge field are linked with a transverse size Ay z ~ ^AAy\. Since A » A r , it 
follows that the characteristic transverse size related with the radiation field 
contribution is much smaller than that related with the space-charge field 
contribution, -\/AA w /2 «: Ay z . By inspection of Eq. (67) and Eq. (68) one can 
see that the value of |r*± - r" ± | is limited by a ± , because of the presence of 
the exponential functions under the integration sign. Therefore, assuming 
constant total charge of the beam, when the electron beam transverse size 
o ± increases beyond ^IAA W the radiation contribution is suppressed with 
respect to the space-charge one. 

Summing up, when condition (2) is valid together with (1), i.e. o]_ » AA W , 
we may neglect the real part of the impedance Z R and approximate the total 
impedance with 



Z = -i 



2coz 



J dr ± J dr" ± p* (r\)p (r 7> ' x )x 



r' 



Ay z 



(90) 



This means that, in the limit o]_ » AA W , the only field to be accounted 
for when calculating impedance (and wake), is the effective longitudinal 
space-charge field. 



5.5 Discussion 



Results obtained in Section 5.4, namely Eq. (90), mean that in the limit o\ » 
AA W radiation is suppressed, so that the beam can be considered as non- 
radiating, and only space-charge impedance is present. Such impedance 
amounts to the free-space impedance, where y is consistently substituted 
with y z . Eq. (90) gives the correct impedance at position z inside the undu- 
lator, as an asymptotic limit for o]_ » AA W of our general theory. 



Our results are in contrast with |13|] - Authors of [130 first noted, correctly, 
that the presence of a finite transverse dimension of the beam o ± suppresses 
radiation in the far zone. Thus, the real part of the renormalized impedance 
in Eq. (80) (valid in the one-dimensional limit of a pencil beam) can be 
generalized to the case when a finite transverse dimension of the beam 
o ± is present, by multiplying it by an exponentially suppressing factor. 
However, they extended such understanding to the imaginary part as well, 
that is incorrect. This led them to obtain the following expression for the 
renormalized undulator impedance accounting for a finite o±_ (see Eq. (2) of 
reference H3]): 
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Z CT± — Z ren exp 
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4cAy : 
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K 2 z lA_ 
2c Ay 2 n \A r j 







-—k a 2 ' 


)] 


exp 






. c 



(91) 



the explicit expression for Z ren being given by Eq. (80). 



Eq. (91) is not derived, in [131 ], as the asymptotic result from a complete the- 
ory It is the outcome of an analogy with the fact that radiation in the far zone 
is suppressed when a finite transverse electron beam size o ± is considered. 
However, the imaginary part of Z ren results as a combination (a difference, 
actually) of two logarithmic contributions, that can be respectively ascribed 
to space-charge and radiative field (see Section 5.1) and present different 
characteristic transverse scales (y z A and y/AA w ) as a consequence of differ- 
ent field formation-lengths. Z ren is calculated in the limit for o]_ <tc AA W . Only 
in this limit the dependence on o ± in the imaginary part of Z ren is cancelled 
as a result of the combination of the before-mentioned logarithmic contri- 
butions. For finite transverse size o±_ such compensation does not take place 
at all. For example, for the ESASE scheme at LCLS (see Section 7), one has 
o ± = 30 /jm, ylAA w 6 fim and y z A = 500/jm. As a result, o ± is small with 
respect to y z A, but large with respect to ylAA w . The incorrectness of Eq. (91) 
(i.e. Eq. (2) of fll3|1 ) follows from these observations. 



Authors of reference [|l3|l conclude that the impedance in Eq. (91), that is re- 
lated with curved trajectory, is suppressed. Thus, the space-charge induced, 
free-space impedance (Eq. (5) in [13]) is finally considered when calculat- 
ing the energy spread inside the undulator. This result is counterintuitive. 
According to it, when the electron beam does not radiate (a 2 ± » AA W ), 
the presence of the undulator does not influence the impedance, indepen- 
dently of the value of the undulator parameter K. However, its incorrectness 
was not trivial to prove, because it relies on an apparently correct analogy 
between real and imaginary part of the impedance. Only developing a com- 
plete theory the presence of two separate logarithmic dependencies can be 
spotted. Thus, in Eq. (90) we saw, as an asymptotic case of our theory, that 
only space-charge impedance is relevant for o\ » AA W , but we additionally 
demonstrated that y must be consistently substituted with y z . 



Our conclusion is that when K 2 > 1 the presence of the undulator strongly 
influences the longitudinal impedance, whether the beam radiates or not 
(i.e. independently of the transverse size o ± ). 
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6 Analytical expression for the wake function in the steady state case 

for o 2 , » AA W 



As discussed in the previous Sections, our derivations hold under conditions 
(1) and drastically simplify under condition (2), i.e. for o]_ » AA W . In fact, as 
we have seen before, when condition (2) holds transverse contributions to 
impedance and wake are negligible. In the following we will consider the 
case when both (1) and (2) are satisfied, that is the case for ESASE schemes 
at LCLS, as we will see in Section 7. 

We will consider a transverse longitudinal profile, as specified in Eq. (72), 
and a longitudinal bunch profile, specified by Eq. (74). Note that the rms 
bunch length a z is connected to the rms bunch duration o t by o z = f}co t , so 
that, in terms of time and frequency we have 



/(*) 



HON 

V2ro f 



exp 



2a 



i — > 



f(co) = (-e)Nexp 



2^21 



CO 



(92) 



When a]_ » AA W , an expression for the wake can be found by Fourier- 
transforming the impedance given in Eq. (90)0. As already noted, Eq. (90) 
is mathematically identical to the free-space expression where only y has 
been substituted by y z . Here we will present only the final result for the 
wake function. For mathematical details regarding wake calculations, we 
refer the interested reader to a previous work of us | j2ll ]. That paper dealt 
with a different subject, namely wake fields and impedances for electron 
beams accelerated through ultra-high field gradients. However, in [21!], we 
also analyzed the steady-state (z » 2y 2 a z ), free-space case of Gaussian 
transverse and longitudinal distribution for the bearr© We fincCE that the 
antisymmetric part of the wake G A is given by G^(A£) = y z r}z/o x ■ H A (A£,), 
where 



H A (AZ) = 



2yfn 



(A£) 2 



iAii 



7i exp 



[(A£) 2 | 


erfc 




4 




[ 2 J) 



(93) 



11 In the more general case when only conditions (1) hold, we should Fourier- 
transform both Eq. (67) and Eq. (68). 

12 In the limit for z » 1y\o z , the antisymmetric part of the longitudinal wake 
function Ga (always defined, as in Eq. (58), from a source disk to a test disk separated 
of As) is dominant with respect to the symmetric part. In this paper we will only 
analyze this part, that will be used later on to discuss the feasibility of ESASE 
schemes. 

13 See Eq. (21) of reference J21Q, where y has been substituted with y z . 
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Fig. 1. Plot of the universal function H A as a function of A£. 

Here we defined A£ = y z (As)/a_ L , t] = y z o z lo± and z = z/(2yla z ). A plot of 
the universal function H A , that is the symmetric part of the wake in units of 
f z z/o±, as a function of A£ is given in Fig. 1. 

The energy change of a single particle at position s within the bunch due 
to the reactive part of the wake, averaged over transverse coordinates is 

Xoo 
m Ga(As)/(s - As)d(As). An explicit 

expression for AS A /S , with S = ym e c 2 , as a function of £ = y z s/cr ± : 





^(0 = ^2*^ fd(A4)H A (£-A4)exp 



(A^) 2 



(94) 



Note that Eq. (94) is a function of £ but also depends parametrically on r], 
and may be presented as 



AS 



7W"r -w F k'V- (95) 

where we indicated the parametric dependence of t] after the semicolon and 
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-15 



Fig. 2. Plot of F in Eq. (96) as a function of s/o z for different values of r\. 

oo 

m 2 ' 



(96) 



A plot of F is given as a function of s/o z in Fig. 2 for different values of rj. 



7 Application to ESASE schemes 



We can now give a practical example of application of our work. Namely, 
we calculate the impact of longitudinal wake fields in ESASE schemes |12l 
13 , li}] . Here we propose an analysis on a set of parameters referring to the 
LCLS i] setup considered in |14Q . Similar calculations may be performed 
on other parameter sets like those for the European XFEL ||9J]. 



We consider a beam with normalized emittance after the dispersive section 
e n - 1.2 mm mrad (like in Fig. 3 of tl3J]). We take the average betatron 
function in the focusing lattice /?/ = 18 m, and y = 2.8 • 10 4 . This gives a 
transverse beam size o ± = (e n fifly) 111 30 jim. The longitudinal size of 
the bunch is o z = 50 nm. The maximal current is about the Alfven current 
Ia — 1 7 kA; in fact, /peak - 18k A. Finally, the undula tor has a period A lv = 0.03 
m, K = 3.7, and the vacuum chamber dimension is a = 2.5 mm. 



We consider a wavelength A a z = 50 nm. We can neglect the vacuum 
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chamber influence (see Section 2), because y z A = 500 fim, as y z 10 4 , and 
y z A <sc a = 2.5 mm. The overtaking length is 2Ay 2 — 10 m. The saturation 
length is about L s = 50 m. Thus z = 5, according to the definition in Section 
6 and we can use our asymptotic expression. Moreover r\ = y z a z /o ± 16.7. 
From fig. 2 (or from direct calculations) one can see that the maximal value 
assumed by F(s/o z/ if) for r\ = 16.7 is about F max - 6. It follows that the 
energy-chirp peak-to-peak is given by (see Eq. (94) or Eq. (95)) : 



A<5 A , P eak = 2m e c 2 ^zF max - 30 MeV . (97) 

LA 

In contrast to this, estimations in lli}] indicate "a swing in energy of 2.4 
MeV". The reason for this large discrepancy is due to the fact that in reference 



114]], where it is correctly recognized that the "most significant cause for 
concern is the longitudinal space charge forces", the Lorentz factor y is 
incorrectly used in place of y z . In fact, as it is clearly stated in that reference 
: "While this expressiorf^ 4 ! has been derived for beam lines containing only 
drift sections and focusing elements, we apply it without modification to the 
present case where the electron beam is passing through the undulator and 
oscillates almost rigidly with a deviation of less than 1/jm". 

In addition to this, it should be noted that energy chirp is also accumulated in 
the free-space between the dispersive section and the undulator, worsening 
the situation even more. In the LCLS case | |l3|] , the dispersive section is a 
dogleg located about 200 m from the undulator. One should account for 
the energy chirp accumulated in this region too, and sum it to that in Eq. 
(97). To this extent, reference |21|1 can be used. The overtaking length is now 
2Ay 2 80 m, so that z = 2.5 and our asymptotic expression for the wake are 
still valid with some accuracy. Using the same procedure as for the wake 
inside the undulator (but considering y instead of y z ), we obtain an extra 
energy chirp of about Afi^peak - 20 MeV. 

The sum of contributions from the straight section after the dogleg and 
from the undulator amounts to about 50 MeV. Although the energy chirp is 
non-linear, in order to estimate the magnitude of the effect we can use the 



linear energy chirp parameter a defined in [15, 16]. The effect of linear energy 
chirp starts to play a significant role on the FEL gain when a > 1. Intuitively, 
this means that the relative energy change becomes comparable with the 
FEL parameter on the scale of the coherence length. The chirp parameter 
is defined as a = -(ya)p 2 1D )~ l ■ dy/dt, p 1D being the one-dimensional p- 
parameter in FEL theory defined as (see ||29j]): 



14 



Eq. (3) of reference \ l4\ . 
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(98) 



where j is the beam current density and the coupling factor Ajj, for a planar 
undulator, is given by A n = J (Q) - /i(Q), where Q = K 2 /(2 + K 2 ). For ESASE 
schemes at LCLS 7 pea k = 18 kA, and we have pm — 10" 3 . Using an estimated 
peak-to-peak chirp of 50 MeV we obtain a — 1. Thus, the saturation length 
is significantly modified |16| ]. This is a reason of concern, because ESASE 
schemes are based on the assumption that the nominal saturation length of 
about 80 m is shortened to about 50 m, that is only 37.5% less. The effect 
described here is fundamental, in the sense that it cannot be avoided by fine 
tuning of the setup parameters. 

Finally, it should be noted that in this paper we did not account for the 
symmetric part of the wake, related with Transition Undulator Radiation. 
This part constitutes only a correction to our calculations, as the space- 
charge wake accumulates along the longitudinal axis, being proportional to 
z. We did not include it in this article, because the space-charge wake alone 
is enough to raise concern. Numerical estimations presented in this paper 
indicate that effects of energy chirp induces by space-charge longitudinal 
wake pose a serious threat to the operation of ESASE schemes at LCLS. 



8 Conclusions 

In this paper we presented a theory of wake fields in an XFEL system, with 



We worked with specific constraints on parameters, that are fulfilled in XFEL 
setups (see Section 2). Namely, we neglected the influence of the vacuum 
chamber and we assumed that the saturation length is long with respect 
to the overtaking length. Our results are valid for arbitrary values of the 
undulator parameter K and in the long wavelength asymptotic, i.e. A » K r , 
K r being the reduced wavelength of the fundamental harmonic. Note that, 
for any FEL setup, the lasing part of the bunch is always much longer than 
A r so that condition A » A r is very natural. It follows that our results are 
of practical importance not only in relation with ESASE schemes, but for 
any FEL setup. We derived expressions for the steady state impedance, 
that is composed of a radiative and a space-charge part. Radiation field 
and space-charge field are characterized by different formation lengths: 
the undulator period K w and the overtaking length 2Ay 2 , respectively. As a 
result, the steady state radiative part of the impedance can be applied for any 
undulator system (with N w » 1), whereas the steady state space-charge part 
of the impedance can be used only assuming that the saturation length is 
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long with respect to the overtaking length, which limits its practical region 
of applicability. Non-steady state results for the space-charge part of the 
impedance can be obtained applying methods presented in |21|1 and |30|]. 



After having dealt with a generic expression for the steady-state impedance, 
we specialized our theory to the case when the transverse beam size o\ » 
AA W . Major simplifications arise in this case: in particular, space-charge con- 
tributions to impedance and wake dominate with respect to radiative con- 
tributions. In this particular condition, that is practically fulfilled for ESASE 
XFEL setups, we showed that the (antisymmetric) wake can be given in 
terms of an asymptotic expression for the wake generated by a beam in 
uniform motion along the longitudinal axis (see ||21|]), provided that the 
Lorentz factor y is consistently substituted with the average longitudinal 
Lorentz factor y z . Final expressions are presented in the case of a planar 
undulator. However, there are no specific effects related with such choice, 
and our work may be straightforwardly extended to the case of a helical 
undulator as well. 

We applied our theory to calculate the effects of longitudinal wake fields 
on ESASE schemes. Our conclusion is that longitudinal wake fields pose a 
threat to the practical realization of ESASE schemes. This finding is in con- 
trast with estimations in literature, where no important detrimental effect 
is foreseen. The reason for this contrast is an incorrect application, in liter- 
ature, of expressions that are valid for beam lines containing drift sections 
and focusing elements to describe the case of XFEL undulators, where the 
longitudinal Lorentz factor is sensibly different. 
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